Heisenberg Spins on an Elastic Torus Section* 
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Classical Heisenberg spins in the continuum limit (i.e. the nonlinear a-model) are studied on 
an elastic torus section with homogeneous boundary conditions. The corresponding rigid model 
exhibits topological soliton configurations with geometrical frustration due to the torus eccentricity. 
Assuming small and smooth deformations allows to find shapes of the elastic support by relaxing the 
rigidity constraint: an inhomogeneous Lame equation arises. Finally, this leads to a novel geometric 
effect: a global shrinking with swellings. ©1998 Elsevier Science B. V. 
PACS: 75.10.Hk, 75.80.+q, 75.60.Ch, ll.10.Lm 
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Curved magnetic structures abound in nature: either 
made out of magnetic materials or enclosing magnetorhe- 
ological fluids. To explore their magnetoelastic proper- 
ties, it is convenient to treat their surfaces as a continuum 
of classical spins. Keeping this in mind, we investigate 
classical Heisenberg-coupled spins on a deformable torus 
section in the presence of topological spin solitons. 

The continuum limit of the Heisenberg Hamiltonian for 
classical ferromagnets or antiferromagnets for isotropic 
spin-spin coupling is the nonlinear a- model Q-^]. The 
total Hamiltonian for a deformable, magnetoelastically 
coupled manifold is given by Ji = 7i m agn +7~tei + 7~L m -ei, 
where Hmagn, 7~lei and H m - e i represent the magnetic, 
elastic and magnetoelastic energy, respectively. In the 
present paper we will focus on the magnetic part and 
the elastic part only. For the nonlinear cr- model, the 
magnetic energy on a curved surface S, in curvilinear 
coordinates, is given by gj^] 



Hi, 



jJJ y/gdSl <r'h.,. .<),„■■;>, 



(1) 



where J denotes the coupling energy between neighbor- 
ing spins. The order parameter n is the local magneti- 
zation unit vector specified by a point on the sphere S 2 . 
The metric tensors (gij) and (h a p) describe respectively 
the support surface S and the order parameter manifold: 
as customary, dfi represents the surface area and g the 
determinant det(g'y). 

First, let us consider the nonlinear cr-model on a rigid 
torus section. The most natural representation of a torus 
is given in cylindrical coordinates (p,£,z): 



p = R + r cos ip, 



r sinp, 



(2) 



where the rotating radius R and the axial radius r must 
verify < r < R, while the angle tp varies from — 7r to n. 
Nevertheless, for our purposes we will use a more suitable 
representation H 



a sinh b 



a suit] 



cosh b — cos T) ' 



cosh b — cos 77 



(3) 



where the new constant parameters a and b are both real 
and positive, while the new angle 77 varies from — ir to n. 
The relations 

R 



a = y/(R + r)(R-r) and cosh b = — (4) 

r 

allow a simple geometrical interpretation for the new pa- 
rameters: let us call a the geometric radius and b the 
eccentric angle. Conversely, the natural parameters R 
and r satisfy 



R = 



and 



tanh b sinh b 

Furthermore, the transformation yields 

tan \r\ = tanh |6 tan \ip, 



R-r 



R 



■ tan ±(p. 



(5) 

(6a) 
(6b) 



One can easily check that the metric on a torus in peri- 
polar coordinates (£, 77) is given by 



9 



(cosh 6 — COS77) 2 
therefore g^ n = g n ^ = and we have 
1 



sinh 2 b d£<g)d£ + dn®dn] , (7) 



9^9- 



sinh b ' 



9 m V9 = sinh 6. 



(8) 
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From now on, we restrict ourselves only to a section of 
the torus: the angle £ about the z-axis will vary from 
— A£ to A£ where < A£ < tt. 

As usual, the local magnetization n is described by its 
polar coordinates (0, $), then the metric on the Heisen- 
berg sphere is given by 



h = de<g)d9 + sin 2 6 d$(g)d$. 



(9) 



Assuming homogeneous boundary conditions at both 
ends (0 = [w] as £ — ► ±A£) allows to map each bound- 
ary of the section to a point: thus we compactify our 
torus section into the sphere S 2 . Consequently, the map- 
ping of our support to the order parameter manifold is 
classified by the homotopy group n 2 (S 2 ) which is iso- 
morphic to Z: spin configurations may be classified ac- 
cording to their homotopy class |l],|9| . 

Henceforth, without loss of generality, only toroidal 
symmetric configurations (d v Q = — 0) will be con- 
sidered. Thus the magnetic Hamiltonian (Q) becomes 



j J d:/ 



S- 

sinh b 



sinh b sin 2 6 <i> 2 



(10) 



where a subscript stands for differentiation. Rescaling 
the angle about z-axis in equation (10) gives: 



+ AC +7T 



n magn = j Jd( fdr) [e 2 + sin 2 e $ 2 ] , (11) 

-AC -7T 

where C = sinh& £ and A£ = sinh b A£. The Euler- 
Lagrange equations corresponding to (O) are: 



® m = 0, 

e fC = $ 2 sine cos e. 



From (12a), it follows that 

*r, = q v q v e Z. 



(12a) 
(12b) 



(13a) 



Substituting this into (12b) and rescaling again the ro- 
tating angle, we get the sine-Gordon (sg) equation 



sin 8 cos 0, 



(13b) 



where g = q ^, o r g = q n sinh b £. 

Equation (13t) may be integrated once to yield 



6 2 = sin 2 9 + m m £ [0, +oo) . (14) 

The limit case m = corresponds to the self-duality 
equation . Performing the change of variable sin O = 
dn(u | 1 + m) where dn is a Jacobi elliptic function [ p^[ , 
the differential equation becomes u 2 = 1. Let us denote 



by a(- | m) the increasing solution of (13b ) specified by 



the parameter m and subject to the boundary condition 
a(0 | m) = |. Readily, we get: 



sina(g | m) — An{g \ 1 + m) . 
Simple calculations lead to the cosine variant: 
cos a(g | m) 



(15) 



-VT+to sn(g | 1 + m) , (16a) 
= - sn^Vl + m | 1/1+m) . (16b) 

Thus, the general SG solution, in natural coordinates, is 

9(g) = ea(g | m) + c e = ±1. (17) 

Further the function a(- \ m) satisfies 

a(g + 2qK a \ m) — a(g \ rh) + qw q G Z, (18) 

where the quasi quarter-period K a is related to the com- 
plete elliptic integral of the first kind K by JlCII 



K a {m) = K(l + m) 
1 



Ki 



- rn 



1 



1 + m 



(19a) 
(19b) 



Using the solutions of equation (13b), the <7£7r-soliton 
configuration consistent with the boundary conditions 
can be obtained easily. Up to an irrelevant additive mul- 
tiple of 7r, we have 

6(£) = ea(q v sinhfe £ | rh) - 5 cvcn . q ^, (20a) 

where the parameter rh is given by 

m = K- l \&A£ sinh6^ . (20b) 

The magnetic energy E mog „ of the above configuration 
( pp| ) may be compared with the corresponding topolog- 
ical minimum energy E which does not depend on the 
geometry of the support manifold. Performing the Bogo- 
mol'nyi's decomposition M yields 



E = 8irJ\Q\, 



(21) 



where the topological charge (i.e. the winding number) 
Q equals to q^q v - Let £ m agn denote the ratio E mas „/E. 
A straightforward calculation shows that £ m agn depends 
on the parameter rh only; we have 



£magn = E(l + m) + \m K(l + m) , 



E 



1 



1 + m 



K 



1 + TTi V 1 + m 



1 



(22a) 
, (22b) 



which increases strictly from 1 to oo with respect to m; E 
is the complete elliptic integral of the second kind Jl(j] . To 
avoid any unnecessary complication, we will consider the 



quantity q v A^/q^ fixed. Therewith, according to (20b), 
£magn is a decreasing function of the eccentric angle b. 
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Clearly, the minimum energy E is only reached when b 
tends to infinity, i.e., when our support becomes an in- 
finite rigid cylinder. This is in agreement with results 
obtained for the rigid cylinder ^,0. In other words, the 
non-satisfaction of the Bogomol'nyi's inequality is due 
to the geometry of the support manifold, hence the ex- 
pression: geometrical frustration. Therefore, here, the 
geometric frustration is induced by introducing a second 
non-vanishing local curvature. Furthermore, notice that 
only the eccentric angle, which measures the balance be- 
tween the rotating radius R and the axial radius r, does 
tune the spread of the soliton. 

Next, let us relax the rigidity constraint and consider 
the nontrivial spin configuration on an elastic torus sec- 
tion. Accordingly, the soliton will try to minimize its 
magnetic energy E mag „ to the minimum energy E by de- 
forming the elastic support 

Since the eccentric angle b appears as the relevant geo- 
metric parameter and the rotating angle £ as the relevant 
curvilinear coordinate, we relax b with respect to £ and 
write 



Before deriving the Euler-Lagrange equation for the 
total Hamiltonian 7i — H ma gn + 7~tei , we calculate the 
magnetic energy associated with the nontrivial spin con- 
figuration (po]). Expanding to second order in A and 
the relations ( 24 ) enables to rewrite (11) as follows 



(1 + ±A 2 ) [9 2 + sin 2 9 $ 2 ] 



T-C magn = J J d( Jdrj 

-AC -is 

-(coth&o A + ±A 2 ) [8 2 - sin 2 9 $ 2 ] 

A 2 e 2 



sinh 6q 



(27) 



On the other hand, according to (|l4|) and (|l5|), the spin 
configuration (EG) verifies 



sin 2 6 $ 2 = mq 2 , 



) q 2 [1 - sn^.C | 1 + m)\ 



b(o = b + HO, 



(23) 



e 2 = (i 



Inserting these relations into (|27| ) and simplifying, we ob- 
tain 



where bo represents the spontaneous eccentric angle and 
the function A describes local deformations. The metric 
on the deformable manifold in this case remains orthog- 
onal, and (||) reads 



1 



sinh b 



, g*yg^wah 2 b + K\. (24) 



As the problem is quasi-one-dimensional, the magne- 
toelastic energy Ti m - e i merely renormalizes the spin cou- 
pling energy J in the magnetic energy TC magn |L2| . There- 
fore, we add to the nonlinear er-model Hamiltonian (jl]) 
only the elastic energy which is essentially stored in the 
bending of the deformable support P^-p^] : 



n el = \k c JJ ygdfl (H — Hq) 



(25) 



Here the constant k c denotes the bending rigidity, H rep- 
resents the mean curvature p7| , and the spontaneous 
mean curvature Hq tends to bias the mean curvature for 
recovering the spontaneous shape. Assuming small and 
smooth deformations and expanding to second order in 
A, A ? and A 55 lead to Q 



n mag n= 2nJtfi J d( (2 + m) - m [coth6 A + ±A 2 ( ] 



-AC 

-(1 + rh) [2 + coth 2 b A 2 ] sn^C | 1 + m) 
1 + m 



||2. 



sinh bo 



A 1 



(28) 



The Euler-Lagrange equation for the problem takes the 
following form 



A ee + [(1 + m)A - mB sn 2 (g | m)] A 

where we have set 

m = 1 + m, 
1 



.4 



nig? 



1 



1 



2 + to V sinh b 



_C2_ 

K(? 2 



B = 2 



J = 



coth bo 
coth bo 



q 2 v VT 



m 



to j, (29a) 

(29b) 
(29c) 

(29d) 
(29e) 



where 



C2 = 



+AC 

H e i = ±nk c c 2 J dCA 2 



-AC 



1 + sinh bo cosh bo (3 + sinh 2 &o) 



sinh b 



Here we have introduced the relative coupling energy 
(26a) K = J/k c and the variable g = q v ( (see (13b)). The 
linear inhomogeneous second-order differential equation 
( p9a| ) is related to the well-known homogeneous (Jaco- 
bian) Lame's equation which occurs in several physical 
contexts ]l8],[l9]]. We have found no direct treatment of 
(26b) fl29a| ) in the literature. However, an approach based on 
the derivation of the Lame functions ITglfill allows to find 
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a particular solution denoted by £(■ | m; A, B,j). There- 
fore, if small and smooth deformations are assumed, a 
suitable deformation function A is given by 



A(£) = £( 9 , sinh ba£\l + m;A,B,j) 



(30) 



The geometric frustration becomes evident in consid- 
ering the Euler-Lagrange equation (29a): A = (i.e. the 
rigid torus section) is not a solution. Now let us turn 
our attention to the deformation mechanism. As we have 
seen, the magnetic energy E magn of the configuration j20| ) 
is a decreasing function of the eccentric angle b: there- 
fore, for sufficiently small and smooth deformations, the 
soliton tries to increase the eccentric angle b. From a 
physical point of view and according to ([|), the soliton 
tends to collapse the torus section (R — > a + and r — > + 
as b — ► oo). On the other hand, the elastic Hamiltonian 
( p5| ) tries to maintain the spontaneous shape of the de- 
formable support: the bending rigidity tends to "attract" 
the eccentric angle b to the spontaneous eccentric angle 
£>o- The underlying physics makes a good sense: rigid- 
ity prevents the support from a collapse ... To sum up, 
any alteration of the eccentric angle b increases the elas- 
tic energy whereas any increasing (resp. decreasing) of 
the eccentric angle b decreases (resp. increases) the mag- 
netic energy — at least for small and smooth deforma- 
tions. Consequently, the competition between the mag- 
netic energy and the elastic energy induces an augmen- 
tation of the eccentric angle b. Further more, according 
to ( |TT| ) the soliton energy is essentially localized in the 
spread zone: the alteration is less important in the spin 
flip region. 

To understand the geometric meaning of the frustra- 
tion release, let us introduce the outer radius R — R + r 
and the inner radius R = R — r. Their geometrical in- 
terpretation is transparent. Now, noting that RR = a 2 , 
we define the relative dilatation A by 

where Ro and R are related to the spontaneous shape. 
The transformation relationships (|5|) allows to write a 
readable relation between the relative dilatation A and 
the eccentric angle b: 



A = 



tanh ±6 
tanh i6 



(32) 



Clearly, the augmentation of the eccentric angle b leads 
to a global shrinking whereas a local swelling arises where 
the spins twist, as shown in Fig. |l|. 



In conclusion, we have shown that the Euler-Lagrange 
equation for the nonlinear cr-model on a torus sec- 
tion with homogeneous boundary conditions is the sine- 
Gordon equation: the spread of the topological soliton 
configurations is tuned by the torus eccentricity. Besides, 
the model presents a geometrical frustration induced by 
the torus eccentricity. As announced, the frustration re- 
lease is shown to lead to a novel geometric effect: the 
torus section is globally shrunk and a swelling appears in 
the region of the soliton. The full torus problem will be 
considered in a subsequent paper flT8| . 

We acknowledge fruitful discussions with J. Mourad. 
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FIG. 1. The relative dilatation A(£) as denned in (^lj) corresponding to the deformation function ( J30| ) associated with the 
7r-Clifford torus section (2A£ = 7r and sinhfeo = 1) in presence of a 7r-soliton (|o|) versus the rotating angle £ for different values 
of the relative coupling energy k = J /k c : the case k = corresponds to the rigid torus section. 
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